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Abstract 

We studied the influence of the amplitude fluctuations of a non-Fermi su- 
perconductor on the energy spectrum of the 2D Anderson non-Fermi system. 
The classical fluctuations give a temperature dependence in the pseudogap 
induced in the fermionic excitations. 
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INTRODUCTION 



The microscopic description of the superconducting state in cuprate materials is a very 
difficult problem because at the present time is generally accepted that in the normal state the 
elementary excitations are not described by the Fermi liquid theory. However, using the BCS- 
like pairing model the Gorkov equations have been applied to describe the superconducting 
state in the hypothesis that the normal state is a non-Fermi liquid described by the Anderson 
mo The superconducting state properties have been discussed by different authorsi~i 
and even if these descriptions are phenomenological, they can be a valid starting point for a 
microscopic model. Recent experimental data (ARPES) showed that these materials present 
even more remarkable deviations from the Fermi liquid behavior due to the occurrence of 
the pseudogap at the Fermi surface. 

The occurrence of the pseudogap has been explained using different concepts as: the spin 
fluctuations!, preformed pairs0, SO (5) symmetryS, spin-charge separation^, the fluctua- 
tions of the order parameter induced pseudogapEl. 

In this paper we start with a non-Fermi liquid description of the superconducting state 
(See Ref. Q-^) and consider the interaction between the order parameter fluctuations and 
the electrons (Section 2 and Section 3). This problem has been studied by Abrahams et 
al.@, Marceljai and Schmid@ for BCS superconductors and the theory explained the tun- 
neling experiments on films, by the modification of the density of states by a pseudogap 
which appears at a temperature higher than the BCS critical temperature. Using such an 
approximation we will calculate (Section 4) the pseudogap due to the electron-fluctuation 
interaction and in the simple mode-mode approximation the temperature dependence of it 
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will be obtained. 

Finally (Section 5) we compare our results with the other theoretical models for the 
cuprate superconductors. 



THE MODEL 

The non-Fermi behavior of the normal state for the cuprate superconductors proposed by 
Anderson0 was developed by different authorsi~i in order to describe the superconducting 
state in the framework of the BCS theory. In the normal state the electrons are described 
by the Green's function 

where uj c is a cutoff energy and < a < 1. 

In the following we consider that the superconducting state appears due an attractive 
interaction and is described by the BCS like order parameter which can be calculated 
from the Gorkov equations. The fluctuations of this parameter can interact with the electrons 
and the fermionic spectrum of the elementary excitations changes. Such an effect has been 
studied in the BCS superconductors by different authorsEUll and it was showed that this 
interaction gives a contribution to the density of states for T > T c which explained the 
behavior of the tunneling measurements. 

For a superconductor described by the Gorkov like equations with the normal state 
described by Eq. (|IJ) the propagator of the fluctuations has the expression: 



D(q,iu n ) = -— — - , (2) 
V 1 + II(q, iuj n ) 



3 



where V is the attractive interaction between the electrons and Il(q, iu n ) is the polarization 
operator defined as 

n(q, iu n ) = t [ 7^^ g (p> ^z)G(q - P> «*>„ - m) (3) 

where G(p, iuji) is the Green's function related to electrons, which in terms of a Dyson 
equation has the following form 

G~ x {p, iuji) = G 1 (p, iuji) - E(p, iui) (4) 

where the self energy is given by 

E(p, iui) = -TJ2J j0j2 D (^ iuJ n)G(q - p, iu n - iui) (5) 

Eqs. (0-0) have to be solved self consistent, but this cannot be done analytically. However, 
in the mode-coupling approximation it can be done and we can calculate the new energy of 
the electronic excitations. 

MODE-COUPLING APPROXIMATION 

In this approximation we consider first that G(\t,iu n ) ~ Go(k, iu n ) and from Eq. (|3|) we 
define the polarization 

f dk 

U (q,iuj m ) = J -^-^S(k,q,iu m ) (6) 

where 

We performed the analytical calculation of n (q, ioj m ) given by Eq. @ (See Appendix) and 
from Eq. ([J) the propagator for the order parameter fluctuations has been obtained as 



D - 1 (q,^ n )=iV(0)A(a)|C(a) 
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where the critical temperature T c has been obtained!!! 
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as 



t: 



2t\ 



C(a) 



A(a)N(0)V 



and the constants from Eqs. @ and (^) are 
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C(2-2a) (-) 2a 



5(3 -2a, a -2) fuo D \ 2a - 2 ( T 
4(2a - 2) V ~u~J \u~ c 



B(x,y) = T(x)T(y)/T(x + y) and T(x) are the Euler' functions and ((x) is the Riemann 
function. 

Using a similar form with the one introduced by Schmid the fluctuation propagator will 
be written as 



I>o (q, iuJ n ) = N(0) b(a)r(a) + ia{a)u n + £ 2 (a, T)q z 
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where 



T 



2a 



r(a) = (-) -1 (12) 



a(a) = [l-a)A{a) (13) 



> rp v 2a 

b(a) = A(a)C(a)[—) (14) 



and 



(ia) = ^lz^l N L2L^) Aia) (15) 



In the approximation £ <C 7rT the Green function given by Eq. @ will be approximated 



we 



as G = Gq + Go^Go an d n will be modified by 511 also linear in S. Following Ref. [1 1 
calculated 511 in the "box approximation" as 

5U = 2T 2 £ / ^GSCp, iu n )Gl(-p, -iu> n ) [ j%D( q , u n = 0) (16) 

n J I 271 ") ■/ (27r)^ 

where 

D _1 (q, iu n ) = — + Il(q, iu n ) + 5II(q, z'o; n ) (17) 
In order to calculate 511 we introduce 

B ° = l (^ G 0(P,^n)G2(-p, _ iUn) (lg) 

where iV(0) = m/27r. If we use for the electronic Green function Eq. ([!]) we obtained 

£(1/2, 3/2 - go) [2 3 " 4 " - 1] C(3 - Aa) 
Bo{T) = n 2^ (vrT)2-4« (19) 

If we introduce f(a) = r(a) + 5II/iV(0) the fluctuation propagator given by Eq. (|17D will be 



D _1 (q, iu n ) = b(a)f(a) + ia(a)u n + £ 2 {a, T)q 2 (20) 



where 



T{a)-T{a)= T / 21 

iV(0) J (27r) 2 r(at) + £ (a, T)g 2 

If we perform this integral taking the upper limit qu = 1 / ^ (ct, T^) from Eq. ( pT| ) we get 

-/ n /x 5 (T)T , l + f(a) 

f (a) - r(a) = °\ J In 22 

For realistic parameters (T c = lOOfT, u; c = 200K) the difference r(a;) — t(q) becomes im- 
portant only near a critical value of a defined by £(cn c ) = 0. In the BCS limit (a = 0) this 
parameter is small and this behavior can be associated with the occurrence of the preformed 
pairs in the domain T c < T < T*, controlled by a. This behavior is in fact due to the 
occurrence of a pseudogap in the electronic excitations. 

ELECTRONIC SELF-ENERGY 

The self-energy due to the interaction between electrons and fluctuations is given by 
Eq. (H) where D(q,iu n ) is given by Eq. (PUD- First we calculate the summation over the 
Matsubara frequencies ui n 

s = TJ2D(i"nMiu n - = — — — — -~r^ (23) 

transforming this sum in a contour integral which has a pole at fi(q) = — (bf + £ 2 q 2 )/a 
and a cut line from + iuj\ to oo in the upper semiplane. From Eq. flT3|) we can see that 
a(a) = —\a(a)\ and in fact O(q) = (bf + C, 2 q 2 )/\a\. Performing this integral we obtain 
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sin W l - a)) r ; /(t) 



_ a;-" sin[7r(l-g)] y°° 

iV(O) 7T 4 (6f- laKt + iwO + W^-ek) 1 " 

where is the Bose-Einstein function and f(x) is the Fermi-Dirac function and e k = 
k 2 /2m — E'p. The integral from the second contribution in Eq. ( |2~3] ) will be performed using 
the expansion f(t) = J2m=o(~ l) m ex P [ — /^( m + 1)*] anc ^ the ^ as ^ term becomes 

"/3(m+l)(fi(q)-e k )' 



_ ~ (-1)- ( £k + f]( q ))^+i 
1 ^ |a| [/3(m + l)]°/2 6XP 
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m=0 

x r(a)WL a/2 ,a/2-i/2 [(3(m + l)(fi(q) + e k )] (25) 

where the Whittaker function Wa,^(<2) will be approximated as W\ tfl = e~ z / 2 z x . This results 
give for Eq. (pl|) the expression 

qg» n(n(q)) 

JV( ) (_^_ £k _fi( q ))i-« 

sin [tt(1 - oQ] ^ (-ir ek + ^ + n(q) 
In the limit k = kp the second term denoted by becomes 

s 2 = ^; ^[(i-^]^H+^+ev r(a)[1 _ 2l _ Q]c(a) (27) 



and if T — > T c , cj; — > and g — ► this term can be neglected. This approximation is 
in fact equivalent with the physical picture proposed by Vilk and TremblayS in which the 
occurrence of the pseudogap is given by the interaction between the electrons and the classical 
fluctuations. Indeed, in this regime the first term of Eq. fl26|) can be written as 

S = ^_n(n(q))G?(k, -M + 0(q)) (28) 
and the electronic self-energy becomes 



E(p, u + zO) = -A 2 pg G(k, -m) (29) 
where we considered 3> fi(q) and 

= TTT^n I ?S;n(n(q)) (so) 



will be approximated as 

a 2 m - T f M qdq ,< v ) 

P9K } 4vr|a|iV(0) h (f + ^V)/|a| 



where qu is the wave number cutoff. From Eq. fl5Ip we calculate the temperature dependence 
of A P g(T) as 

A » (T ) = 4^ to ( 1 + T^) (32) 
DISCUSSIONS 

We showed that a temperature dependent pseudogap appears in a non-Fermi supercon- 
ductor due to the interaction between electrons and the fluctuations of the order parameter 
amplitude. 

The mode-mode coupling, valid in the weak coupling approximation, can give relevant 
results, even for the intermediate coupling studied by the Levin grou p0 usine; the resonant 
scattering model. The method, recently applied by Norman et a 10, RanderiaS can be applied 
for the spin-fluctuation model proposed by Chubukovi in order to study the temperature 



dependence of the pseudogap. In Ref. [18] and [19] the filling in of the pseudogaps due to 
the incresment of the temperature is given by the broadening in the self-energy and is 
proportional to T — T c . A similar broadening effect, proportional to f(a) was obtained in 



our model and this can be seen very easy from Eq. fl27D if in the electronic Green function 
we take the limit q = 0. 

Recently such a model, for a Fermi liquid superconductor has been studied by Kristoffel 
and Ordfl and their temperature dependence is different from our result. However, we men- 
tion that according to their model these authors have to obtain a result similar to the result 
given inHH. The difference is given by the method of performing the integral over q which is 
not correct in@. 

Recently, Preotsi et al0 generalized the method given inlll taking into consideration the 
anysotropy in the dynamic susceptibility due to the interplane pairing. 
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APPENDIX 



The polarization Il(q, iu m ) for a 2D non-Fermi liquid is defined as 

r dk 

n (q, iuJ m ) =TJ2 ~pr\2 G o(^ iv n )G (q_ - k, iu m - iu) n ) 

= T y[^ ^ (33) 



which can be written as 



f dk 

n (q,i^ m ) = J j^r^S(k,q,iu m ) (34) 
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where 



-2a 



S(k, q, iuj m ) = (-l) 1_t *T Y - — — -2 : — - (35) 

In order to perform the summation in Eq. ( |35|) we transform he summation in a contour 
integral 

S(k,c l ,z) = -f c ^-n(z)F(z) (36) 
where n(z) is the Fermi function and F(z) is given by 

U-a, .-2a 

F W = 7 ^r 1 — -^r- ( 37 ) 

and the contour C is taken as (— oo, iaj n + £ q _k) U(^k? oo). The integral in Eq. (|36|) has been 
evaluated as 



dz „f~\i?t~\ - 1 / P £q " k ^„^ , \ 2it«;- 2Q sin[7r(l -a)] 



/ (12 , , n , , 1 /— £ q-k 

<p nmFm = < / dxnix + iu) 

Jc2ni K ' K J 2ni \ J-oo K 



(x - e v ) 1 - a {x + iu n + eq-k) 1 



, , , 2iu- 2a sin vr 1 - a) 1 , N 

dxn(x)- ^ ^— — \ 38 

In order to perform the integral over x we express the dominators from (^) as 

(x - e k ) a ~\x + £q _ k + iu n ) a - 1 = (x - 5k) - 1 (x + ek)^ 1 

— (a — l)(t)_p5 cos 9 — ic<j n )(x — £k) a_1 (a ; + ^k) a2 

+ (( " ~ (v F qcos6) 2 (x - e k ) a ~\x + e k ) a " 3 (39) 



and take for the Fermi function the expansion 



oo 

n ( x ) = J2 (-l) m exp [-/3(m + l)x] (40) 

m=0 



Using now the integrals 
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dx(x + (3)- v (x - uY~ l = (u + Pf- v B(v - fj,, n) 



dx{x + f3f u -\x - u) 2 "- 1 exp [-fix] 



{u + /3) 



// 



v+p 



exp 



(P - u)fi 



x T(2p)W u ^ P)U+p ^ 1/2 (ufi + Pn) 



we calculated S(k, q, iu m ) as 



S(k, q, LO m ) 



u c 2a sin [7r(l — a)] 



7T 



X 



{-(2e^f a - 1 B(l-2a ) a) 



+ 2T(a) g (-1)^ ^ a _ 1/2M(m + 1)] 



m=0 



(m+ l)]«-i/2- 



+ (a- 1) (v F q cos 6 -uj m ) [-(2e- k ) 2a ^B(2 - 2a,a - 1) 
1 10 ^ £ (~l) m r^^^ ^-3/ 2 [^/3(m + 1)] 



m=0 



(m+ l)]°-3/2- 



+ 



(a — \) 2 {v F qcos6)' A 



(2e k ) 2a ~ 3 B(3-2a,a-2) 



P( m + !)] 7 
2r(a-l) ~ (2 £k ) a ~ 3/2 ^ ' 



7T 



m=0 



(m + 1)] Q - 3 / 2 ' 



Eq. fl34|) will be written as 



n(q, z^ m ) = 2N(0) / — / deS(e, q, iw. 

JO 27T JO 



and using the relation 



x»K v (ax)dx = 2^ 1 a~'"- 1 r ( = ' C ' ' ) T 



1 + // + zA /! + // — 1/ 



we obtained 



n (q,w) 



2iV sin [vr(l - a)] 



7T 
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+ 



[ 2a 

iu(l - a) r o2a ~ 2 



2a 



UJr 



2 2a ~ 2 B{2 - 2a, a - 1) /u^X 2 "- 1 
2a - 1 W / 



r(a - l)r(g - 1/2) 1 - 2 2 - 2 % ,„ / T 

+ f= ~o . 2ct C(2a - 1) 



'7T 



2 2 - 



OJ, 
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+ 
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2 2Q " 3 J B(3-2a,a-2) /cuzA 2 "- 2 
2a - 2 V ~uj 



(v F q) 2 (l - a)' 
Au 2 



1 - 2 



3-2q 



2 3-2a 



C(2-2a) (-) 



2a-2H 



(45) 



Using now the Thouless criterion 



1 + VReU(q = 0, iu m = 0) = 



(46) 



we calculate the critical temperature 



2a 



B(l-2a,a) ( uj d\ 
2aT 2 (a)(l-2 l - 2a )C(2a) \~J 



2a 



7T 



N{0)VT 2 {a)2 2a {l - 2!- 2 «)C(2a) sin [tt(1 - a)] 



(47) 



which is identical to Eq. (R) if we introduce A = N(0)V. 
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